REPORT No. 793 



EXPERIMENTS ON DRAG OF REVOLVING DISKS, CYLINDERS, AND STREAMLINE RODS AT 

HIGH SPEEDS 

By Thbodobh Thbodorsbn and Arthur Rbqier 



SUMMARY 

An experimental investigation concerned primarily with the 
extension of test data on the drag of revolving disks, cylinders, 
and streamline rods to high Mach numbers and Reynolds num- 
bers is presented. A Mach number of 2.7 was reached for 
revolving rods with Freon 113 as the medium. The tests on 
disks extended to a Reynolds number of 7,000,000. Parts of 
the study are devoted to a reexamination of the von Kdrmdn- 
Prandtl logarithmic resistance law and the AcJceret-Taylor 
supersonic drag formula and conditions for their validity. 
The tests confirm, in general, earlier theories and add certain 
new results. A finding of first importance is that the shin 
friction does not depend on the Mach number. Of interest, 
also, are experimental results on revolving rods at very high 
Mach numbers, which show drag curves of the type familiar 
from ballistics. A new result which may have general applic- 
ability is that the effect of surface roughness involves two distinct 
parameters, particle size and particle unit density. The part- 
icle size uniquely determines the Reynolds number at which the 
effect of the roughness first appears, whereas the particle unit 
density determines the behavior of the drag coefficient at higher 
Reynolds numbers. Beyond the critical Reynolds number at 
which the roughness effect appears, the drag coefficient is found 
to be a function of unit density. In the limiting case of particle 
lt saturation,^ or a maximum density of particles, the drag 
coefficient remains constant as the Reynolds number is 
increased. 

THEORETICAL BACKGROUND 
von kArman-prandtl theort for pipes 

Measurements of the value of the skin friction between a 
fluid and a solid constitute one of the means for studying 
the nature of turbulent flow. Most of the pioneer analytical 
work in this field is found in the papers by von K&rmfin (ref- 
erences 1 and 2) and Prandtl (reference 3). The treat- 
ment used in the first part of this section follows the work of 
Prandtl which, in turn, is closely related to the von K&rm&n 
papers. The theory, which concerns the flow in pipes, is 
given in considerable detail as it forms the basis for the suc- 
ceeding discussion on flat plates, cylinders, and disks. The 
theoretical work in this section constitutes mainly an at- 
tempt to analyze and organize earlier work found in many 
scattered articles. Considerable work along such lines has 
already been done by Goldstein, who is responsible for an 
expression for the drag on revolving disks. 

The von K&rm&n-Prandtl theory for flow in the turbulent 
layer is based on the following two assumptions: 



(1) The ratio of the velocity deficiency to the friction 
velocity is a function of geometric parameters only. 

(2) Adjacent to the wall, but beyond the laminar sublayer, 
the slope of the curve representing this ratio is inversely 
proportional to the distance from the wall. The constant 
of proportionality is a universal constant. 

The friction velocity is defined as 

and the corresponding friction length is defined as 

(All symbols used in this paper are defined in appendix A.) 
A reference time may be given as 

T _L_ V _fl 

The geometric conditions for a pipe are given by one param- 
eter, the radius a. A revolving cylinder of infinite length 
represents another single-parameter case, in which the refer- 
ence parameter is the radius of the cylinder. 
The equation of motion can be written in the form 



u , ( y t a\ 



and, by adopting suitably defined mean values with respect 
to time, at a given profile 



u - (y a\ 



Henceforth u will designate such mean velocity. By 
measuring the velocity with respect to a velocity U e in a 
fixed geometrical position c=ka, 



u 



is obtained. About 1930 von K&rm&n showed that for the 
turbulent layer this function is essentially independent of L 
and dependent only on the geometry as indicated in assump- 
tion (1) ; therefore 

u-U c _ 
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This quite remarkable relationship, which has been generally 
confirmed by Nikuradse, Wattendorf , and others (references 
5 to 7), implies a similarity in the turbulent-field pattern 
away from the walls at all Reynolds numbers. The basic 
reason for this similarity remains unknown. 

It follows from assumption (2) that near the wall 

u-U * 



=- log U+ Constant 



where I Ik is the constant of proportionality. (Natural 
logarithm has been used throughout except where otherwise 
indicated.) Since u=U s at y=5, this relation reduces to 

This logarithmic relationship holds to a certain value c of 
the significant parameter a (see fig. 1), where c=ka with k a 
constant. The value of 1— k is only a small fraction, so 
that the point c will be relatively close to the wall. The 
velocity in the center of the pipe is therefore given as the 
sum of three expressions, that is, 



For the laminar sublayer 



s 

- z =a 



and the equation may be rewritten as 

%-«->«+^H>°e|+[Xf)I 



where 
and 



Ci = a log a 



The constant Ci is equal to the nondimensional velocity 
measured on the logarithmic velocity profile when this 
curve is extrapolated to y=L, and the constant & is the 
excess velocity in the center of the pipe as compared with 
that of the logarithmic line extended to y=a. (See fig. 1 .) 
When these constants are combined, the following general 
relation is obtained: 

Umax n _,l 1„„ <*> 

The application of this theory to cases other than circular 
pipes is restricted to geometric configurations given by a 
single parameter. Tt is interesting to observe that both 



Ci and 1/k are universal constants resulting from the second 
assumption — namely, that the flow near a wall is a function 
of the distance from the wall only The second constant 
Ci which gives the excess velocity as compared with the 
logarithmic distribution at a reference point, the location 
of which depends on the geometric dimensions involved, is 
not a universal constant but is dependent on the configuration 
and the choice of reference length. 

The effect of surface roughness may be treated in a similar 
manner. If the roughness parameter e/L is less than a certain 
magnitude, there is obviously no effect at all. This value 



of e/L is found experimentally to be 3.3. 



For I 



L > 3 * 3 ' 

Umax/Ur is shown to be constant, or independent of L, 
except for the so-called unsaturated condition which will be 
defined later. Thus 

^=(7+1 log 3.3^ 

=<7+ilog3.3+ilogf 

or 

%=tf 1+ Ilogf 

The velocity distribution is exactly as if there were a 
laminar layer present of a thickness 8 « 3.5 e or as if the length 

L were ^ e- When ^^3^3 € > the velocity distribution no 

longer changes with an increase in Reynolds number R. It 
seems, therefore, that the distance from the wall of the 
innermost disturbance, or the mean value of the thickness 



I8r 



Universal velocity 
deficiency function. 

a 




Figure 1. — Parameters and functions of the velocity profile by the von Kannan-Prandtl 

theory. 
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of the laminar layer, is of the order of three to four times 
the height of the irregularities or the grain size €. This fact is 
not inconsistent with the physical interpretation. 

The quantity U^Ur is shown to equal J— Further, 



and, therefore, 



_ V Umax 



where B is referred to the maximum velocity and is equal to 
Umaxd/v. The equation 



Umax n \ 1 a 

_=<?+- log z 



may thus be written 



or 



where 



o 3 =- 



V2 



By the similarity hypothesis, the mean velocity in a pipe 
differs from the maximum value by a constant, or 



U m 



u a 



T7T 



where U m is the mean value of the velocity. Prandtl gives 
4.07 for the value of K 2 . (See reference 3, p. 142.) Note 
further that the product R-yfO^ remains the same whether 
B and 0 D refer to the mean or the maximum value of the 
velocity; therefore, 

^=tf-4.07+i log 
and, finally, with B and C D referring to the mean velocity, 

where 



With 0=5.5 and k=0.4, 



0—4.07— log -42 



V2 



This value is not accurately estabhshed, as the various 
authors seem to differ. 

DRAG OF FLAT PLATES 

In order to obtain the drag formula for flat plates, a 
calculation similar to the von Karman-Prandtl treatment 
for pipes may be performed. The velocity deficiency Au is 
given by the relation 

where U Tm is a mean value between 0 and x y the distance 
along the plate. The missing momentum may be written as 



or 



= J 0 u dy -l Ku) dy 



M 

w 

where U is the stream velocity and Si is a significant length 
giving the thickness of the boundary layer. Rewritten, 
this equation becomes 

or, by virtue of the similarity law, 

Since the momentum is given directly as 



the following identity is obtained : 
«. 1 V" 



or 



which gives 



2 



% C& - Co / Gp m 



Using the logari thmi c deficiency relation gives for C& the 
value 1/k, or 2.5, and for C 6 f0 6 the value 2/*, or 5; thus 



h i V : 



Dm 



x 2.5 
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By use of the von Kfirmdn-Prandtl treatment, the stream 
velocity is obtained in essentially the same form as for pipes. 
With small adjustments, therefore, 

__* 3+ -l0g-— 

By use of the expression for Si/s, the following equation is 
obtained: 

^-*.+4.071o gl „j-M* te 

LOCAL VALUES OP DRAG COEFFICIENT FOR FLAT PLATES 

It may be noted that a relation for the local drag coefficient 
on a flat plate may be found in a fashion similar to that used 
later for a disk. Consider a plate of unit width; for the full 
length Z, 

d=c d (± p u*i) 



With the subscripts m and x referring to mean and local 
values, respectively, for the length x, 



= J 0 Dx dx 



or 



Bv 



dC, 



Therefore 
where 



R 



<7Dx=C Dm (n+l) 



Dx 



dQogE) 

BOUNDARY RELATION FOR REVOLVING DISKS 

The moment coefficient is defined as 

M 



G M = 



The moment may also be written 
M=2p J (2wa)u,u& dy 



=We.(^)'tX'^Mi) 

where u T is the variable radial velocity and u t is the tangen- 
tial velocity, from which 



or 



-« Constant 



The drag formula then reads 

A similar result was obtained by Goldstein in reference 4. . 
TESTS AND RESULTS 

Tests on disks, cylinders, and streamline rods were con- 
ducted to determine drag or moment coefficients. For the 
cylinder the two coefficients are equivalent; for the disk and 
the rod it is more convenient to employ the moment coeffi- 
cient, which can be measured directly. In order to extend 
the range of Mach number, several tests were conducted 
with Freon 12 or Freon 113 as the medium. The test 
results obtained are of technical interest because some of 
the data, particularly for the high Mach number range, 
were obtained for the first time. It may be pointed out 
that many of the earlier tests on revolving disks and, in 
particular, on revolving cylinders were conducted on a 
rather smali scale and in a limited range of Reynolds num- 
ber. It may be noted that a considerable range of Reynolds 
number is generally needed in order to confirm with suffi- 
cient reliability a particular theoretical formula. For 
instance, it may be impossible to obtain a measurable 
difference between logarithmic or power formulas if a short 
range of Reynolds number is available. This matter of 
distinguishing between the various types of formulas is of 
theoretical interest. 

EXPERIMENTS ON REVOLVING DISKS 

The moment coefficient is defined as 

a M 

— 1 

This definition corresponds to the one for laminar flow on a 
revolving disk given by von Kirm&n in reference 1 as: 



where 
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The constant a* used by von K&rmfin was 1.84 for one side 
or 3.68 for both sides; this value was later adjusted by 
Cochran (see reference 8, voL I, p. 112) to ai=3.87. If this 
corrected value of a x is inserted, the formula for laminar 
flow reads 

Ck=3.87i2- 1/2 

The turbulent-flow formula as given by von Kdrmfin for 
revolving disks is 

Oir=0.146B- 1/5 

In figure 2 are shown the experimental results for tests of a 
series of revolving disks. The Reynolds number ranged 
from about 1600 to more than 1,000,000. Note that the test 
points lie along the theoretical curves given by the von 
Kdrmdn formulas. The transition from laminar flow is seen 
to occur at jR=310,000. This was the largest value reached 
with the most highly polished disk. 

The thickness of the laminar boundary layer is, according 
to von Kfinndn, 



5=2.58. 




or, which is equivalent, 

i=2.58fl- 1/2 
a 

Using 22^=—^ leads to 



§=1=2.585-^ 

For the transition Reynolds number, 310,000, 

.Rj=2.58V/B 
= 1440 

which is of the same order as the minimum critical value 
obtained for pipes. 

Several tests were conducted for the purpose of investi- 
gating the factors affecting the transition Reynolds number. 
The first observation was that the transition Reynolds num- 
ber could not be increased beyond the value 310,000 no 
matter how highly the surface was polished or whatever 
other precautions were taken. Likewise, it was unexpectedly 
difficult to decrease the transition Reynolds number. The 
application of coarse sand (60 mesh) glued to the surface of 
a disk (1-ft radius) only reduced the transition Reynolds 
number to about 220,000 (fig. 2). The reduction in the 
transition Reynolds number by initial turbulence was also 
studied. A small high-pressure air jet applied near the cen- 
ter of the disk produced the greatest observed reduction 
(fig. 2) and brought the transition to a point near the inter- 
section of the lines representing the drag formulas for lami- 
nar and turbulent flow, which is the absolute minimum. 
Note that the drag in the turbulent region is quite appre- 
ciably increased by surface roughness. 




3.6 3A 



4.8 5.0 52 
Log„R 

Figure 2. — Moment coefficient Cm~~^~ — for disks as function of Reynolds number. 
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The values of the moment coefficient given in figure 2 
represent obviously an integrated drag over the disk. An 
expression may be obtained for the local drag coefficient C D x 
as a function of local Reynolds number as follows: 



=2 £ Cds G pwV ) (27rr2)rfr 




dC 2 
i 



•dO : 



By substituting 



2?r eta 47r 



or 
and 

where 
If 

then 



C ^"27r G3f L^aogii!) + 2j 



log (7Dx=log Ck+log -^r 

T r_ d(log C M ) 
d(log22) 

C M =cB M 

n _ 5+2n r 



4tt 

By use of the expression for log 6j>*, some of the data of 
figure 2 are plotted in figure 3. Although the general picture 
does not change much, the abrupt nature of the transition 
becomes apparent. 

An illustration of the boundary-layer profiles for various 
radii or Reynolds numbers is given in figure 4, in which 
curves of equal velocity ujosr are also plotted. Note that 
the thickness of the boundary layer in the laminar region is 
essentially constant. The transition value of R, 310,000, 
is shown approximately by the line marked "Approx. 
transition" in figure 4. The nominal laminar boundary- 
layer thickness consistently appears to be somewhat in excess 
of that given by von K&rm&n in reference 1. There appears 
to be some discrepancy from the theoretical velocity distri- 
bution which is shown for the laminar boundary layer as 



-L6 




Figure 3.— Local drag coefficient for disks obtained by differentiating the moment-coefficient carve. 
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obtained from work by Cochran. (See reference 8, vol. I, 
p. 112.) It is recognized that the experimental error in this 
case is of considerable magnitude. The turbulent boundary 
layer shows almost perfect agreement with the logarithmic 
curve, which is plotted for one profile in figure 4. 

It may be remarked here that a series of hot-wire tests 
were run to study fluctuations in the boundary layer with the 
following results: 

(1) No disturbances were noted in the laminar region 

(2) A pure tone of a frequency of about 200 cycles per 
second was observed in the transition region 

(3) A random disturbance involving much higher fre- 
quencies was observed in the turbulent region 

In figure 5 the upper range of the Reynolds number has 
been considerably extended. The highest Reynolds number 

reached is 7,000,000. The y-power law holds fairly well in the 

observed range which, however, is too limited to permit a 
distinction between the power law and the logarithmic law 
for the velocity distribution. The main purpose of the tests, 
the results of which are shown in figure 5, was to investigate 
the effect of the Mach number. The first run taken with air 
as the medium extended to a Reynolds number of about 
2,000,000 and a Mach number of 0.62. By using Freon 12 
as the medium, the range of Reynolds number was extended 
to 7,000,000. At the lowest pressure, the highest value of 
the Mach number reached was 1.69. All the data for Freon 12 
show a slightly higher drag than that given by the von 
Karman formula, apparently because of some systematic 
error. The significant result of this investigation is that the 



drag coefficient is absolutely independent of the Mach 
number. A separate extension of the experiments to a Mach 
number of slightly more than 2 further confirmed this inde- 
pendence of the Mach number. 

EXPERIMENTS ON REVOLVING CYLINDERS 

The experimental results for revolving cylinders are shown 

in figure 6 as a plot of logiod? against logiofi, where R=^-* 

The drag formula for laminar flow on a revolving cylinder 
is obtained from Lamb (reference 9, p. 588) as 



Cd ~r 



where 



D 

M 
Ta 



In this formula S is the surface area and a the radius. In 
this case it is convenient to use C D instead of C M} which was 
used for the revolving disk, because no integration is in- 
volved. The laminar curve is shown in figure 6. The drag 
relation given by 

1= -0.6+4.07 log 10 fiV<k 

ye© 

for the turbulent flow is also shown in figure 6. 



// Disk radius 
405000 R 




O A .8 

Fiqueb 4. — Observed velocity profiles on revolving disks showing the transition red on. 
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<6 48 50 52 

Figure 5.— Moment coefficient for disks as function of Reynolds number for several Tames of Mach number with air and Freon 13 as mediums. Maximum Maeb number, 1,00. 
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The experimental results are replotted in figure 7, where 
is shown as a function of logioR-y/0 D . The relation 
for the turbulent flow 

1 =-0.6+4.07 logMUv 

appears in figure 7 as a straight line. The coefficient 0 D in 
this formula corresponds to a value of 0.4 for von Karman's 
universal constant k. The relation for the laminar region 
4 

£#=22 appears as a curved line near the origin. 

It is noted that the drag coefficient for rough cylinders is 
dependent on the relative grain size e/a } where e is the size 
of the sand and a is the radius of the cylinder (see fig. 8), 
and that for each grain size the drag coefficient remains 
constant and independent of the Reynolds number beyond 
a certain minimum or critical value, which lies on the line 
for turbulent flow. In regard to the magnitude of the drag 
coefficient as a function of relative grain size for particle 
"saturation" of the surface, it may be remarked that the 
value of 5 is a measure of the thickness of the sublayer or, 
what amounts to the same thing, a measure of the minimum 
grain size of the turbulence. It is therefore to be expected 
that the surface roughness will become effective at the 
Reynolds number for which €„, the critical value of c, be- 
comes less than the grain size e. Inversely, it may be seen 
that, if the Reynolds number becomes smaller than this 
critical value, the grain size of the turbulence is too large to 
be affected by the surface roughness. With € greater than 
€„, which is 3.3Z, the following relation is approximately 
true for the drag coefficient beyond the critical Reynolds 




is 24 32 

Figure 7. — Drag parameter for smooth cylinders as function of logtoA Vcp* 
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Figure 8.— K fleet of surface roughness, or grain siie, on the drag coefficient. Saturation density of particles. 
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number for surface roughness of saturation density: 
-4=== -0.6+4.07 logio 3.3 V2- 

=2.12+4.07 log 10 - 

In figure 9 the experimental points are shown to satisfy this 
theoretical relation with sufficient accuracy. 

Tests were made to determine the effect of the density of 
spacing of grains of a given size, and the results are presented 
in figure 10. Such tests were made with a certain unit 
grain size but with the surface density in grains per square 
inch varied between 90 and 2200. The grain size used'cor- 

responds to the size ^=0.03, also used for the preceding ex- 

(if 

perimental results shown in figure 8. It is verified that the 
critical Reynolds number depends on the grain size only, 
and it is further shown that the slope of the drag curve 
beyond the critical Reynolds number is a function of the 
density. A saturation condition evidently always exists, 
in which the drag coefficient remains approximately con- 
stant and equal to the critical value. 

EXPERIMENTS ON STREAMLINE RODS 

In figure 1 1 results are given for certain more or less stream- 
line bodies, each tested in two or more different mediums. 
The tests were obtained by using actual propellers of 12-inch 




*4 .8 12 1.6 2D 24 28 32 3j8 4D 42 48 52 56 GO 

Figure 10. — Effect of varying density of surface roughness on the drag coefficient. 
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diameter, which are designated propellers B and C. Pro- 
peller B had a section of double synunetry with a circular-arc 
contour line. Propeller C was obtained by reducing the 
chord of propeller B by removal of about one-fifth of the 
chord near one extremity to obtain a blunt-nose airfoiL By 
running propeller C backwards an airfoil with a blunt trailing 
edge could also be studied. The coefficient used in figures 
11, 12, and 13 is the standard torque coefficient for 
propellers 



For the symmetrical airfoil B, a value of the Mach num- 
ber 1 of about one was reached in air, the range was extended 
to 1.6 in Freon 12, and the characteristic decrease in the drag 
coefficient was finally reached in Freon 113. A- considerable 
decrease in drag coefficient was noted at the largest Mach 
number, 2.7, which to the knowledge of the authors is the 
highest Mach number reached except for a few cases of 
projectiles. 
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Figure 11, — Torque coefficient Cfr— ^ipi M faction of Mach number for propellers B and C. 
i Note that the Mach numbers used in figures U, 12, and 13 are based on the tip radius. 
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Figure 12.— Torque coefficient Oq" J^j^ as function of Mach number for propeller E. 
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Figure 13.— Torque coefficient Cfr— as function of Mach number for propeller D. 
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The blunt-nose airfoil section C showed approximately the 
same low-speed resistance as the symmetrical sharp-nose 
section B but had a maximum torque coefficient very much 
in excess of that of section B. The test extended only to 
near the peak of the torque curve with Freon 12 as the 
medium. By reversing the direction of motion of propeller 
C to obtain a blunt rear, the expected large increase in drag 
at low Mach numbers was observed. The appreciable 
difference in Reynolds number for air and Freon 12 is 
apparent from the difference in drag coefficients in the range 
below a Mach number of unity. For higher Mach numbers, 
the drag coefficient of the section with the blunt rear lies 
between the drag coefficients of the doubly streamline section 
and the blunt-nose type; the streamline leading edge is 
approximately twice as effective as the streamline trailing 
edge, a result in general agreement with earlier observations. 
It should be noted, however, that the lowest drag is obtained 
with both leading and trailing edges streamlined. 

The effect of the Reynolds number is also shown in figure 
12, which gives the results of tests to study how the scale 
effect is superimposed on the Mach number effect. It should 
be noted again that the Reynolds number effect appears only 
for a Mach number below unity. A wide variation in the 
Reynolds number shows no consistent measurable effect on 
the drag for a Mach number greater than unity. Similar 
data for a small angle of attack, instead of zero angle of 
attack as used in the preceding discussion, were used in one 
case, for which results are given in figure 13. 

The four propellers referred to in figures 11 to 13 are shown 
in a photograph (fig. 14) and the dimensions of the propellers 
are given in table I. 

TABLE I 

DIMENSIONS OF PROPELLERS OR REVOLVING RODS FOR 
TESTS AT HIGH MACH NUMBERS 



All propellers have a straight taper In chord and thickness. The tips are rounded as shown 

in fig. 14.] 



Propeller 
designation 


Airfoil section 


Pitch 
(deg) 


At 50 percent radius 


At 92 percent radius 


Chord 
(In.) 


Thickness 
On.) 


Chord 
(In-) 


Thickness 
(in.) 


B 


Circular arc ~ 


0 


L75 


a 31 


L07 


a 14 


0 


Blunt nose 


0 


L30 


.35 


.82 


.15 


D 


Circular arc 


•2.5 


L44 


.18 


L03 


.11 


E 


Circular arc 


0 


.S3 


.13 


.52 


-07 



• Propeller D was twisted so that approximately the outer half of the blade had an angle 
of attack. 



It is of some interest to interject a superficial analysis of 
the results presented herein, in view of Ackeret's formula 
as given by Taylor (reference 10). For the local section 
Ackeret gives the drag coefficient as 

where the bar indicates the mean value. For zero angle of 
attack and a symmetric section with j9 1 2 =)8 2 2 J this relation 
becomes 



For a circular-arc section, 0 2 =- jS Mai 2 , where p nax is the 

maximum angle. This angle is, in turn, approximately 
equal to twice the thickness ratio t, which is the total thick- 
ness divided by the chord. For circular-arc sod ions, 
therefore, 

Figure 15 shows G D plotted against Mach number for 
different values of i. At M— 1.0, the curves tend erroneously 
to infinity. This effect follows from a simplifying assumption 
used in the derivation of Ackeret's formula. 




Figure 14.— Propellers B, 0, D, and B. 



.24 




O A .8 12 1.6 2.0 24 2.8 

Mach numher) M 



Figure 15. — Theoretical curves of the drag coefficient Co against Mach number for various 
thickness ratios for circular-arc airfoils by Ackeret's formula. 
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By using the general form/ (M) instead of the Mach num- 
ber function -5— 1, the drag coefficient may be written 

c D =^ mm 

The torque coefficient is known experimentally to be a 
function of the Mach number, or ljx u where Xi is the fraction 
of radius at which the Mach number is unity; thus, the fol- 
lowing integral relation is obtained: 

There are several ways of handling this relation. The non- 
dimensional chord c and the thickness t may be taken to 
represent a preferred section at approximately 80 percent 
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Fiqxjhe 16.— Values of drag function f(M) as function of Mach number from analysis of 
experimental moment curves for propeller B in figure 11. 

of the radius. By assuming an initial drag coefficient 0 D 
any desired accuracy may be obtained by iteration methods. 

The function f(M) shown in figure 16 has been obtained 
for propeller B by such a process based on the experimental 
data given in figure 11. Note that the drag coefficient 
approaches the value given by the Ackeret formula for large 
values of M, for which /(M) approaches (M*— 1)~ 1/2 . Note 
further that the maximum value of the drag coefficient 
occurs at M=1.2 withjf(M) almost exactly equal to unity. 
It is, of course, not to be concluded that the function f(M) 
has general validity; the function is given here for propeller B 
for the purpose of comparing the data with the Ackeret 
theory, 

CONCLUDING REMARKS 

Experimental results on the drag of revolving disks have 
been presented, which substantiate to a remarkable degree 
drag formulas based on the von Kdrmdn-Prandtl theory of 
skin friction. The range of the investigation was extended 
to a Mach number of 1.69, which is beyond the range of any 



earlier test, and to a Reynolds number of 7,000,000. It was 
established that the skin friction is independent of the Mach 
number up to this value and appears to be a function of the 
Reynolds number only. 

The drag at supersonic speeds was studied with revolving 
rods or propeller sections. Mach numbers as high as 2.7 
were attained in the tests. The drag at supersonic speeds is 
a function of the Mach number only, as it appears to be essen- 
tially independent of both the Reynolds number and the 
nature of the medium. The characteristic peak in the drag 
curve observed for projectiles was obtained. For thin 
streamline bodies, this peak appears at Mach numbers only 
slightly beyond unity; in fact, it appears at a Mach number 
of about 1.2. Systematic tests were conducted on stream- 
line bodies with combinations of sharp and blunt leading 
and trailing edges for the purpose of obtaining the relative 
merits of such features. It was found that the increase in 
the peak value of the drag coefficient resulting from a blunt 
nose is about twice that resulting from a blunt trailing edge, 
when both drag coefficients are compared with the drag 
coefficient of a section with streamline leading and trailing 
edges, which has the lowest value. 

Significant results were obtained on revolving free cylin- 
ders for which references to earlier tests seem to be lacking. 
It was found that, at very low Reynolds numbers, the drag 
asymptotically approaches the laminar drag of the classical 
theory whereas, at higher Reynolds numbers, the drag is 
found to conform to a logarithmic formula of the von 
K&xm&a type. There is no distinct transition from laminar 
to turbulent flow, as is found in pipes and on revolving disks. 
The flow is essentially turbulent down to the smallest 
Reynolds numbers. 

The effect of initial turbulence was particularly studied in 
connection with tests of revolving disks. It was found that 
the transition Reynolds number was very slightly affected. 
The critical Reynolds number at which the roughness effect 
appears depends on particle size only and is not a function 
of particle density. Beyond this value of the Reynolds num- 
ber, the drag coefficient is constant only when the surface 
is "saturated," that is, when the density of the individual 
particles attains a maximum value. For a roughness of less 
than this particle density, the drag coefficient decreases with 
Reynolds number. 

It is interesting further to note the persistence of the 
logarithmic relationship. When 1/V(7d is plotted as a func- 
tion of log R-yfU^ (where G D is the drag coefficient and B is 
the Reynolds number), the lines representing turbulent flow 
are invariably straight. A rather critical demonstration of 
the logarithmic velocity pattern near the surface is thus 
shown. The range investigated is of considerable extent. 
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APPENDIX A 



SYMBOLS 



TO 
P 

u 

v maI 
u m 

v c 

u 

An 

u r 
u t 

03 
d 

*i 
L 
I 

T 
t 



v 
M 

r 
a 



friction velocity ( 



y 

1 

' K 



shear per unit area at surface 
mass of air per unit volume 
mean friction velocity (from 0 to x) 
stream velocity for flat plates 

rnnvimnTn velocity 

mean velocity (in pipes) 

reference velocity (at a given fraction of 

radius or of other reference dimension) 
velocity at 5 

absolute variable velocity of fluid in 

boundary layer 
velocity deficiency, stream velocity minus 

local velocity for flat plates 
radial velocity for disks 
tangential velocity for disks 
angular velocity, radians 
thickness of laminar sublayer 
boundary-layer thickness 
friction length (v/U T ) 
total length of plate 
reference time (LfU T ) 

time; also, thickness ratio for propeller 
/Thickness of airfoiPv 

sectlc H cm ) 

coefficient of kinematic viscosity 

coefficient of viscosity 

variable radius of pipe, disk, or propeller 

radius of pipe, cylinder, or disk; also, 
velocity of sound in fluid 

distance from leading edge of flat plate in 
direction of flow; also, fraction of pro- 
peller radius ( z= ]g where R denotes 
radius of propeller tip^ 

fraction of propeller radius at which Mach 

number is unity 
distance normal to surface 

nondimensional profile constant for turbu- 
lent flow near walls 



fraction of reference dimension 



also, nondimensional chord of airfoil 

/ChordA 

VRadius/ 



Cdz 

D 

D x 

e 

o u 

M 

R 
R» 

R x 



R* 

Re 
V 

2 

S 

o Q 
Q 

N 
n 

Pi, ft 
fima 

(k 



c%, Ct, . . . 

Ki, jKi, 
k 



angle of attack of airfoil; also, profile 
constant (<5/Z) 

total-drag coefficient (Many authors use /, 
7, or X/4 instead of Cd for pipes.) 

mean drag coefficient (from 0 to a;) 

local drag coefficient 

drag; also, propeller diameter 

drag of plate (from 0 to x) 

grain size of roughness 

grain size of critical roughness for particu- 
lar value of drag coefficient 

moment coefficient for revolving disks 

missing momentum; moment for disks; or 
Mach number 

Reynolds number 

Reynolds number based on thickness of 

boundary layer 
Reynolds number based on distance from 

leading edge of flat plate or on local 

radius of disk 
Reynolds number based on pipe diameter 
Reynolds number based on pipe radius 
velocity (Ackeret formula) 

dynamic pressure ^for cylinders, ff =-|pw 2 a 2 ^ 

area of cylinder 

torque coefficient (Qlpn 2 D 5 ) 

torque 

number of blades 

rotational speed, revolutions per second; 

also, coefficient in power law 
angles which upper and lower surfaces of 

airfoil make with center line 
maximum angle which circular-arc section 

makes with center line 
nondimensional velocity measured on log- 
arithmic velocity profile when this curve 

is extrapolated to y=L 
nondimensional excess velocity at y=a 

over that of logarithmic line extended 

to y=a 

constants 
constants 
constant 

constant in equation for moment coeffi- 
cient of revolving disks 
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APPENDIX B 



NUMERICAL VALUES OF POWER REQUIREMENTS FOR REVOLVING DISKS AND CYLINDERS 



A chart is presented (fig. 17) which gives the horsepower 
required to drive a smooth disk in standard air (760 mm 
and 15° C, p=0.00238 slugs/cu ft and u=0.000159 ft'/sec). 
Lines of constant horsepower ranging in value from 0.01 to 
1000 are plotted with disk rotational speed (in rpm) as 
abscissa and disk diameter (in ft) as ordinate. The dashed 
line in figure 17 represents a Keynolds number of about 
400,000, which is considered the transition Keynolds 
number. 

The following formulas were used to calculate the power 
for disks operating in the turbulent region: 



co-war**?- 1 '* 



=0.146 



Horsepower==^jj 

0.146 
^550X2 



Inasmuch as the formula for Cu is based on the 1/7 
power for velocity distribution, the calculated values of Cm 




too 



1,000 10.000 
Rotational speed, rpm 



40,000 



Fioube 17.— Power requirement for smooth disks 
843110—00 20 



are too low for high Reynolds numbers. This error may 
become appreciable for the highest power, since the chart 
(fig. 17) covers a range of Reynolds numbers to 60,000,000. 

A chart is also presented (fig. 18) which gives the horse- 
power required to rotate a smooth cylinder of unit length 
(1 ft) in standard air. The following formulas have been 
used in calculating the curves: 



Mc0=Cx£SO6) 



la< 



CO 



rr Mo 

Horsepower = 

~ 550 
where, for smooth cylinders, 

-0.6+4.07 log 10 RJC D 

and, for rough cylinders, 6>€ er , 

vfe =2 - 12+4 - 071og »r 




too 



WOO 10,000 
Rotational speed, rpm 



^aooo 



Figure 18-— Power requirement for smooth cylinders (Mt length). 
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APPENDIX C 



COLLECTED SKIN-FRICTION FORMULAS 
FLAT PLATES (ONE SIDE) 



SYMBOLS 



The following symbols are used iu the formulas for flat 
plates collected herein: 

0 D total drag coefficient 

0 Dx local drag coefficient at point x 

x distance from leading edge of flat plate in direction of 
flow 

I length of flat plate in direction of flow 

R Reynolds number based on I 

R x Reynolds number based on x 

LAMINAR FLOW 

The formula for total drag coefficient 
<7 i >=1.328fi- 1/2 

is based on the simplified hydrodynamic equations developed 
by Prandtl in 1904. (See reference 2, p. 2.) The constant, 
which was calculated by Blasius in 1908 as 1.327, was cal- 
culated by T6pfer in 1912 as 1.328. (See reference 3, p. 89.) 
The formula for local drag coefficient is 

<7z>*=0.664B x - 1/2 

Von K&rmfin, Schoenherr, and others have indicated that, if 
the total drag coefficient is 

<7x>==Constant R* 

the local drag coefficient is given as 

C D£ = (n+l)0 D 

This relation is derived in the section entitled **Local Values 
of Drag Coefficient for Flat Plates" in this paper. All for- 
mulas given in this appendix for the local drag on flat plates 
are in conformity with this derivation. 



TURBULENT FLOW — SMOOTH SURFACE 



The formulas 
and 



(7z>=0.074B- 1/6 

<7x> x =0.059i5! x - 1/5 

were first calculated by von K&rmfin in 1920. (See references 

1 and 2.) Based on results from pipes and on the y-power 

law for velocity distribution, they are consequently valid 
in. the lower Reynolds number range, ii<10,000,000. 
382 



Some writers use the following formulas of the same 
type, which are fairly accurate to a Reynolds number of 
500,000,000: 

C r D =0.030£" 1 ' 7 

Cz> x =0.026B x - 1/7 

Of more general validity are the so-called logarithmic drag 
formulas of the type 

=4.15log 10 fitfz> 



The form of this relation was determined by von K Arm An 
with constants adjusted to conform with data by Schoenherr 
and' others. (See reference 2, p. 12.) In the present paper 
a different form has been developed, which is in somewhat 
stricter theoretical conformity with the physical relations 
involved: 

1 RO 

veT 4 - 07 logl ° 1-3.MVCS 

Prandtl has developed an explicit expression which gives 
essentially the same results as the logarithmic formulas. It 
is 

C r z>=0.455aog, 0 B)- 2J8 

(See reference 3, p. 153.) The local drag coefficient has also 
been given by von K&rm&n in a logarithmic form with the 
constants adjusted to fit the experiments of Kemp, which 
included measurements on small movable plates inserted on 
a long pontoon. This formula is 



= 1.7+4.15 log 10 fi x <7i>* 



(See reference 2, p. 12.) 

TURBULENT FLOW— ROUGH SURFACE 

Schlichting (see reference 8, p. 382) gives the two following 
formulas for the total and the local drag coefficients for rough 
flat plates, respectively: 



C D =(l.89+1.62log IO 0 
Oz*=(2.87+1.58 logiof)" 



•2.5 



Von K&rm&n (reference 2, p. 18) gives for the local drag 
coefficient for rough surfaces a formula of the logarithmic 
type 

-^==5.8+4.15 log 10 f 



DRAG OF REVOLVING DISKS, CYLINDERS, AND STREAMLINE RODS AT HIGH SPEEDS 

R x Reynolds number at radius x 

LAJVUNAR FLOW 
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PIPES 

SYMBOLS 



The symbol R d used in this section refers to the Reynolds 
number based on the pipe diameter and the mean flow 
velocity, and the symbol R a refers to the Reynolds number 
based on pipe radius. Some writers use / or y instead of 
C D , used herein, and others use X where X=4Cz>. 

LAMINAR FLOW 

For laminar flow in pipes the formula for drag coefficient is 

This formula is attributed to Poiseuille and Wiedeman. 
(See reference 3, p. 38, and reference 8, p. 298.) 

TURBULENT FLOW— SMOOTH SURFACE 

The formula for drag coefficient for turbulent flow in 
smooth pipes is 

O D =0.079Rr 1/4 

This formula is based on the experimental workof Blasius 
(see reference 3, p. 136), for which the Reynolds number 
range was rather limited. Later work by Nikuradse (ref- 
erence 5) extended the range of Reynolds number to a 
much higher value. The following formula of the type 
developed by von Kfirmdn fits the data better: 



-1== -0.40+4.00 log 10 BdV^> 

(See reference 8, p. 338.) In the present paper a formula of 
this type with different constants is developed: 

-4^=0.40+4.07 log 10 R a ^fC~ D 

TURBULENT FLOW— ROUGH SURFACE 

For turbulent flow in rough pipes 

-^=3.46+4.00 logm* 

Tho experimental work in deriving this formula was done 
by Nikuradse. (See reference 8, p. 380, and reference 6.) 

REVOLVING DISKS 

SYMBOLS 

The following symbols are used in the formulas for revolv- 
ing disks: 

Ofs moment coefficient 

Onx local drag coefficient at radius x 



For laminar flow 
and 



Cu =3.877?- 1 2 
O^—Rr 1 " 



This formula for local drag coefficient is derived from the 
relation 

For the development of this relation and for references, see 
the section entitled "Experiments on Revolving Disks" in 
this paper. 

TURBULENT FLOW 



For turbulent flow 



C r a/ =0.146/t- 1 ' 5 



and 



<7,*=0.053fi J - 1/6 



The formula for the local drag coefficient 0 Dx is derived 
from the equation for the moment coefficient G M in the same 
way as for the case of laminar flow. The local drag coefficient 
in logarithmic form may be given as 

-^==-2.05+4.07 log M RtJUd 

The constant —2.05 has been adjusted to fit the data of 
figure 3. 

REVOLVING CYLINDERS 
For laminar flow 



4 
R 



For turbulent flow on smooth cylinders 

^==-0.6+4.07 log 10 iJV^ 

For turbulent flow on rough cylinders 

-^=2.1+4.0 log 10 <* 

The development of these formulas and the references are 
given in the section entitled "Experiments on Revolving 
Cylinders " 
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